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\£) . Abstract 

Vacuum polarization in QED in a background gravitational field induces interactions which 

(N 

effectively modify the classical picture of light rays as the null geodesies of spacetime. After a short 
introduction on the main aspects of the quantum gravitational optics, as a nontrivial example, we 



study this effect in the background of NUT space characterizing the spacetime of a spherical mass 



o 

\ endowed with a gravitomagnetic monopole charge, the so called NUT factor. 
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I. INTRODUCTION 



It is not unexpected that the effect of QED vacuum polaraization in a background grav- 
itational field, though tiny, should modify the photon propagation in that field. Vacuum 
polarization is an effect in which the photon exists as a virtual e + e~ pair for a short time. 
This virtual transition assigns photons with an effective size of 0(X C ), where A c is the Comp- 
ton wavelength of the electron . Obviously the photon propagation will be affected by 
the gravitational field if the scale of the spacetime curvature L is comparable to A c . The 



influence of the vacuum polarization on the propagation of a bundle of rays has 



Deen studied 



a. 



Ignoring 



recently through the perturbative deformation of the Raychaudhuri equation 
the vacuum polarization, the equivalence principle leads to the photon propagation at the 
speed of light on the spacetime null geodesies. But when the pair production is taken into 
account, the equivalence principle is violated in such a way that the photon trajectories are 
modified and superluminal (photon) propagation, without necessarily breaking down the 
causality [2] , becomes a possibility. 

A. Quantum gravitational optics 

Considering the vacuum polarization, the first order correction to the electromagnetic 
action 

1 /v 



W Q = ~J ^x^F^ (1) 
is given by 

W x = ^ j d*x{-g)*{aRFp V F> a ' + bR^F^F» + cR^F^F^ + dD^D v F v a ) (2) 
in which 

la .13 a la ,1a 

d=- — - (3) 
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are perturbative coefficients of 0(a), m e is the electron mass and a is the fine structure 

constant. Expression ([2j) is called the Drummond-Hathrell action [lj]. The first three terms 

represent the influence of curvature while the last one exists even in the flat space-time an 

shows the off-mass-shell effects in the vacuum polarization [lj. The effective equation of 

motion can be obtained by varying W = Wo + W\ with respect to A u as follows, 

D^ v - — {2bR^D»F° v + 4ci?/ ar D >1 F' JT ] = (4) 



in which only terms of first order in a and (X c /L) 2 are kept. Now we use the geometric optics 
approximation in which the electromagnetic field is written as a slowly-varying amplitude 
and a rapidly-varying phase through a small parameter e, i.e, (A^ + ieB^ + ■ ■ ■ )exp(z^), 
where $ is a real scalar field. In geometric optics, the wave vector is defined as the gradient 
of the phase, i.e, k^ = d^d and the polarization vector is introduced through = Aa^, 
where A is the amplitude and is the po 
gives the light cone condition at (9(e _1 ) as 



arization vector. The geometric optics method 

a 



2b v 8c v 

k -R^yk^k -\- —R^i/fjfk^k a a = 0, (5) 



„2 

m~ m z e 



in which a^a^ = — 1. Since this effective equation of motion is homogeneous and quadratic 
in k^ it can be written as 

Q^k^K = 0. (6) 

In other words at this order of approximation there is no dispersion and consequently the 
phase and group velocities are equal and given by 

v P h = |r (7) 

Now there is a possibility that the modified null cone lies outside the original light cone 
leading to superluminal photon velocities. The photon trajectories can be deduced from the 
modified equation of motion as 
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ds 2 aT ds ds ml 



(bR aT — AcR apT \a p a x 



dx a dx 7 



apT ds ds 



0. (8) 



In the classical general relativity the photons move along spacetime geodesies which are 
the integral curve of fc M = ^p. This is no longer the case when quantum corrections are 
taken into account due to the modification produced by the third term in eq. (jSJ). For 
Schwarzschild space-time, it is shown that the paths depend on the polarization of the 
photon and are given by classical trajectories plus a correction term which depends on the 
black hole mass, the electron mass and the photon energy. 



II. QUANTUM CORRECTIONS IN NUT SPACE 



The NUT spacetime line element is given by 



fl: 



ds 2 = e~ 2u [dt - 2q(l + cos 6)d(f)] 2 - (1 - ^)"V^r 2 - r 2 d6 2 - r 2 sin 2 ^0 2 (9) 
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where 



-2u 



i--( 9 2 + mv^; 



(10) 



This is a stationary spherically symmetric solution of the vacuum Einstein field equations 
with two parameters, the mass M and the NUT factor q. Obviously when q vanishes, the 
metric reduces to that of the Schwarzschild spacetime. It is more convenient to introduce a 
local orthonormal frame with the following basis 1-forms 

e° = e ° t dt + e° 

e 1 = e 1 r dr 

e 2 = e 2 e d6 (11) 

e 3 = e 3 A6 



where the vierbeins are 



A 



-2qe~ u (l + cost 



e% = l/Vl-^/rV 
e% = r 
e 3 , = r sin 6 



(12) 



Calculating the connection one forms it could be seen that there are nine non-vanishing 
components of Riemann tensor out of which the following three are independent 



R 



0101 



A -R0202 = B R, 



0312 



C 



(13) 



where 



A = ( 1 - 3 ] e~ 2u I v" - 2v' 2 + 



q V 



r 3 (l-2l 



C = qr -\ 1 - ^-e~ 2u (1 + rv' 



Introducing the following notation 



(14) 



U2, = W-Wa U = 0,1,2,3, 
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the Riemann tensor can be written compactly as: 

R abcd = 2 A(U m U m - U 23 U 23 ) + 2B(U 02 U 02 + U 03 U 03 - U 12 U 12 - U 13 U 13 ) (15) 
+ C(U 03 U 12 + U 12 U 03 - U 02 U 13 - U 13 U 02 - 2U m U 23 - 2U 23 U m ). 

To keep the expression clear we have suppressed the lower indices a, b, c, d in all terms. For 
a Ricci flat spacetime, the light cone condition, eq. (jSJ), is given by: 

k 2 a b + eR abcd k a k c a d = 0. (16) 

where e = — J^-. This is a set of three simultaneous linear equations for independent com- 
ponents of the polarization a b . It is convenient to introduce the three independent linear 
combinations of momentum components as 

h = h a UH m b = k a U° a 2 n b = k a U°J (17) 

together with three dependent combinations 

(18) 



Pb = 


k a U a l = 


1 (it 
^(k m b 


- k%) 


Qb = 


k a Ul 3 = 




- k 3 l b ) 


r b = 


k a U 23 = 


^(k 2 n b - 


- k 3 m b 



The vectors Z, m and n are independent and orthogonal to k a . To solve photon propagation 
equation, we rewrite eq. ( fl6|) as a set of equations for the independent polarization compo- 
nents a.l, a.m and a.n by substituting for the Riemann tensor from eq. (TIB]) . Doing so we 
arrive at; 

k 2 a.l + 2eA (l 2 l.a) + 2eB [m.lm.a + n.ln.a — p.lp.a — q.lq.a — r.lr.a) 
+ eC (n.lp.a + p.ln.a — m.lq.a — q.lm.a — 2l 2 r.a — 2r.ll.a) = 
k 2 a.m + 2eA (l.ml.a) + 2eB (m 2 m.a + n.mn.a — p.mp.a — q.mq.a — r.mr.a) 
+ eC (n.mp.a + p.mn.a — m 2 q.a — q.mm.a — 2l.mr.a — 2r.ml.a) = 
k 2 a.n + 2eA (l.nl.a) + 2eB {m.nm.a + n 2 n.a — p.np.a — q.nq.a — r.nr.a) 

+ eC (n 2 p.a + p.nn.a — m.nq.a — q.nm.a — 2l.nr.a — 2r.nl. a) = (19) 

Although difficult, in principle it is possible to solve the above system of coupled equations 
in general. We select special trajectories which illustrate the most important features. 



A. Radial component of motion 



For the radial component of motion in which the photon momentum components satisfy 
k 2 = k 3 = eq. f[T9~l) could be written in the following matrix form; 

( k 2 + 2eA{-k°k° + k 1 k 1 ) \ //„\ 

k 2 + 2eB(-k°k° + k^k 1 ) 

y k 2 + 2eB{~k°k° + k 1 k 1 ) J 



La 
m.a 
y n.a j 



(20) 

One can solve for the eigenvectors of the above matrix as the polarization vectors satisfying 
(I16p . The k 2 value for each polarization is obtained through the vanishing of the determinant 
of (|20|) . Since e is a small parameter [2J we keep only the terms up to 0(e) and therefore 
all the second terms in the diagonal components are ignored so that the only degenerate 
eigenvalue is A: 2 = 0. This is a result expected from the spherical symmetry of NUT space 
and it means there are no corrections to the radial component of motion and in consequence, 
the photon velocity is given by | p- 1 = 1 for all the polarizations. 



B. Orbital component of motion 

Now consider the orbital component of motion in which the photon polarization is such 
that k 1 = k 2 = 0. Diagonalizing the corresponding matrix we find as eigenvalues Ao = k 2 
corresponding to the unphysical polarization vector, a&, proportional to nj, and 



A+ = k 2 ± e 



9C 2 (k°k 3 ) 2 + (A - Bf (k°k° + k 3 k 3 )' 



(21) 



corresponding to the transverse polarization eigenvectors. From this the velocity shift is 
found to be 



5v = ±--=-e 
2 r 6 



1 + \w (f ) 2 ^ ~ 2 ™ 2 > + Hi (fl'^f + lw ^ 2 + 8M2 ) 



9 1 fq\ 



9 1 /<r 



6 / q 



24 fq- 



(22) 



The pure mass term, = ^-^2^ is factored out to see the contribution of the NUT 



factor in the velocity shift as compared to that of mass. 
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III. THE TRAJECTORIES 



To obtain the new trajectories the modification term in equation (jSj) should be calcu- 
lated along with the classical approximation to 0(a°). Although these equations are too 
complicated to be useful in general, employing the space-time symmetries will help us to 
obtain useful information about the quantum modifications. In this respect we note that the 
third term in (jSJ) is a derivative of a scalar and so the covariant derivative can be replaced 
by a partial derivative. On the other hand, as mentioned before, this scalar is polarization 
dependent but since it is calculated along the classical paths it possesses all the classical 
symmetries. In space-times with spherical symmetry such as NUT, Schwartzschild and 
Reissner-Nordstrom, it depends only on the radial coordinate r and so the only modified 
component is the r-component. Employing the (1+3) threading formulation of space-time 
decomposition f| to NUT space it could be written in the following general form; 



ds z 



e' 2v (dt - A a dx a f - laP dx a dx p . 



In the classical case varying x a , we find 



d 



1 d 



5x a \ -f [e~ 2v (i - ApiP) A a + lap x?\ + (i _ ApxPf - 



(23) 



0. (24) 



dr i v~ -t>~ /-« ■ 2 dx c 

Here the projected 3-dimensional metric 7^/3 can be written in the form involving the unit 
Cartesian vector f 



<y a pdx a dx> 3 = e 2X dr 2 + r 2 (dr 



(25) 



where e 2A = (1 — 1 e 2u and dr the variation of f is an arbitrary small vector perpendicular 
to f. Now making variations with respect to f and r and inserting the modification term 



along r, we obtain 



6] 



d 



dr 



, df\ 



df 



— r — + 2k— x qr = 



dr 



dr 



and 



K 2 e 2u _ .2 e 2A _ L 2 r -2 + J_J_ / ^ R a Md^_ 

mi or \ dr dr 



where L — r x 4- and k = e 2u (i — A a dx a ^j . Crossing both sides of eq. (|26il by 



we end up with 



d f 2 df\ dL d 
rx— r— = — = —2— (nqr) 
dr V dr J dt dr y H ' 



L + 2nqr = J = const. 



(26) 
(27) 

(28) 
(29) 



i.e. the cone underlying the classical null trajectories is not modified by the quantum correc- 
tions. In classical limit, to see the geometry of the trajectory, we introduce the curvilinear 
angle around the cone's surface. Then r 2 = L, so eq. (jUJ) can be integrated by quadrature 

Lr " dr (30) 



^ K 2 e -2(X-u) _ L 2 r -2 e -2X 

In general this integral can not be performed explicitly for A and v of NUT space, however, 
it has a C (r) solution, which solved for r gives r c (0). The corresponding quantum modified 
trajectory is given by 

Lr~ 2 dr 



-2(A- V ) _ L 2 r -2 e -2\ + e ~2\f ( r> a ) 



/(r,a) = i-UcR^a^ d 4-\. (31) 



dr \ dr dr 

with a solution in the form of = C (r) + <f>i (r, a) which can be solved to give r = 
r c (0) + r i (0?o). Trajectory modifications induced by vacuum polarization are along the 
r-component for the two physical polarizations, equal in magnitude but opposite in sign 
{sj]. The two modified trajectories, corresponding to two transverse polarizations, on the 
r — <f> plane, are expected to be shifted in opposite directions with respect to the classical 
path such that the difference from the classic path varies with the value of the r coordinate 
at each point on the path. We recall that the trajectory splitting happens on the classi- 
cal cone and reduce to the space-time null geodesies when a — > 0. The similar situation 
happens in Schwartzchild and Reiner-Nordestrom cases. For Schwartzchild space-time, the 

Q 

corresponding shift for circular trajectories has been demonstrated in [7J. 
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